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to  a  collision,  this  leads  to  a  more  difficult  problem,  not  discussed  in  this  paper.  However,  see 
[Sh]  for  a  treatment  of  the  two-dimensional  version  of  this  more  general  problem  by  a  method 
using  generalized  Voronoi  diagrams. 

In  the  problem  studied  here  we  wish  to  determine  whether  any  pair  of  spheres  in  the  given 
collection  5  intersect.  Note  that  if  instead  of  spheres  we  consider  rectangular  parallelipipeds 
(which  we  will  call  boxes  for  short)  whose  faces  are  parallel  to  the  standard  coordinate  planes, 
then  the  problem  of  detecting  intersection  between  some  pair  of  these  bo.xes  can  be  solved  in  time 
0(n  log-«)  by  using  known  2-D  dynamic  range  query  techniques  (cf.  [Ed]). 

To  handle  the  more  difficult  case  of  spheres,  we  need  the  following  geometric  lemma. 
Leniina  1:  Let  Si,  S2  be  two  intersecting  spheres  having  centers  Cj,  co  and  radii  r,,  r^  respec- 
tively, and  suppose  that  r,  >  rj.  Let  AT  be  a  cube  inscribed  in  »i  and  having  faces  parallel  to  the 
coordinate  planes.  Then  if  the  angle  8  between  CiCo  and  the  line  /  parallel  to  the  2-a.\b  and  pass- 
ing through  Co  is  <oc,  where  sina  =  — -pr,  then  /  intersects  K. 

Proof:  Define  d„^a  ^°  ^'^  t^^  smallest  possible  angle  that  can  be  formed  between  UjUo  and  /  such 
that  I  lies  either  e.xterior  to  A'  or  on  the  boundary  of  K .   Then 

.                     rf(«^         r^/v/3-          r^           1 
sin  flnaln  =   "77 T  -^   -^   ~:; =  ~  =  Sin  a 

Hence  if  ^  <  a  we  abo  must  have  d  <  9^^^,  from  which  it  easily  follows  that  /  intersects  K. 

Q.E.D. 

The  algorithm  we  propose  is  roughly  as  follows.    Let  "I*  be  a  finite  set  of  orientations  having 

the  property  that  for  each  orientation  $  there  exists  another  orientation  ^6*1'  such  that  the 

angle  between  0  and  0  is  <a.    Let  A/  be  the  size  of  'J'.  For  each  Oq^  ^  v/e  apply  the  following 

space-sweeping  technique.  Assume  without  loss  of  generality  that  ^o  points  in  the  direction  of  the 

positive  z-a-\is.    Sweep  a  plane  P  parallel  to  the  x,y  plane  from  :  =  -oo  to  j  =  +  co.  As  P  is 

swept  upwards  we  maintain  a  data  structure  T  (to  be  described  below)  containing  some  of  those 


spheres  in  S  whose  centers  lie  below  P.  Each  time  P  passes  through  the  center  u  of  a  sphere 
5  6  5,  we  perform  the  following  steps: 

(a)  Let  r  be  the  radius  of  s.  Perform  a  range  query  which  enumerates  ail  spheres  a  in  7  whose 
radius  p  is  not  larger  than  r  and  whose  center  projects  into  a  point  of  P  which  belongs  to  the 
square  of  side  2r/v/3  centered  at  u.  Check  each  such  sphere  a  for  intersection  with  a.  If  an  inter- 
section is  found  the  algorithm  reports  this  intersection  and  halts.  Otherwise  delete  the  sphere  a 
from  T. 

(b)  Add  the  sphere  s  to  T . 

It  is  clear  that  the  sweeping  procedure  just  described  performs  n  range  queries  and  n  insertions 
into  T.  Moreover,  each  sphere  a  E  S  can  appear  in  at  most  one  range  query,  since  it  is  deleted 
from  T  immediately  after  any  such  query  Also  the  number  of  sweeps  performed  by  the  algorithm 
is  a  constant  independent  of  n.  We  will  show  that  T  can  be  maintained  in  0(n  log  n)  space,  in 
such  a  way  that  individual  insertions  and  deletions  from  T  require  time  O(log^n),  while  each 
range  query  needed  can  be  carried  out  in  time  O(log^n  +  t),  t  being  the  number  of  spheres  satis- 
fying the  range  condition  of  the  query.  The  preceding  remarks  show  that  use  of  such  a  structure 
T  will  lead  to  an  algorithm  requiring  0(n  log  n)  space  and  0(n  log'n)  time. 

The  data  structure  T  that  we  use  represents  a  collection  8^,  .  .  .  ,3^  of  spheres.  For  each 
such  sphere  s,  let  x,,  y,  be  the  x  and  y  coordinates  of  its  center  and  r,  its  radius.  T  is  simply  a 
balanced  binary  range  tree  with  a  certain  amount  of  auxiliary  data  stored  at  its  nodes.  The 
leaves  of  T  correspond  to  the  spheres  Si,  .  .  .  ,»t,  which  are  sorted  by  the  x  coordinates  of  their 
centers.  Each  internal  node  ^  of  f  represents  an  interval  [^,„„,  ^^,,f,\  of  z-values,  and  to  each  such 
{  we  attach  a  priority  search  tree  Qf  (as  defined  by  McCreight  in  [McC])  which  represents  the 
set  of  all  spheres  v.hose  centers  have  x  coordinates  lying  in  this  range;  each  such  sphere  is 
represented  in  Q^  by  the  pair  \y ,  r]  of  the  y  coordinate  of  its  center  and  its  radius.  As  shown  by 
McCreight,  each  such  priority  search  tree  Q  has  size  0(m)  and  supports  query  operations  of  the 


form  £"(yrir.  i'la,.  r-_„].  which  enumerate  all  pairs  [y.  r]  in  Q  for  which  y-^  <  i'  <  l-rai  -D'^ 
r  <  r_^.  Each  such  qaer>-  requires  time  0[log  m  ~  i).  where  m  is  the  number  of  pairs  in  O. 
and  where  .'  is  the  number  of  pairs  satisfying  the  query  conditions  to  handle  each  query.  Inser- 
tion? and  deletions  of  elements  in  such  a  tree  require  time  0(!of  m)  each  Not-e  that  the  space 
occupied  by  the  range  tree  T  to  whos«  nodes  the^e  priority  search  ;:e«s  are  attached  is 
Oil  log  n). 

It  is  easy  to  apply  the  operations  we  reqnire  to  T.  To  insert  a  sphere  i  (represented  by  a  triple 
[i.  y.  r]  giving  the  r  and  y  coordinate  of  its  center  and  its  radius)  into  T  we  insert  it  into  ea-rh 
of  0(!og  n  I  priority  search  trees  Q:.  for  ^  lying  on  the  path  from  the  root  of  T  to  an  appropriate 
leaf  of  7.  inserting  the  pair  y.  rj  into  each  of  these  Or.  Each  of  these  updating  operation  can  be 
accomplished  in  <?(log  n)  time,  so  that  the  overall  cost  of  inserting  a  sphere  into  T  is  Ollog^r. ). 
Deletions  are  performed  in  a  symmetric  and  sLmilar  fashion.  Each  of  the  queries  that  we  apply  has 
the  form  E.\{z-^.  T;,  y._.  y^.  r-).  and  calls  for  the  enumeration  of  all  spheres  (i.e.  triples  [z.  y.  r') 
in  T  satisfying  'i  <  -  <  -^  ^i  <  y  <  3/;.  :ind  r  <  r,-,.  To  apply  scch  a  query,  we  find 
0(log  ".J  nodes  in  T  representing  disjoint  inten  als  which  collectively  span  the  inter.al  between 
Zi  and  r;.  For  each  such  node  (  we  can  apply  the  qxiery  E'\-j-.  y-,  -;)  to  9.-  in  lime 
O(log  n  +  ir).  thus  enumerating  all  the  '.:  spheres  satisfying  the  conditions  of  the  original  query 
whose  centers  have  ^-coordinates  hing  in  [z.n-  C-.-t'-  Overall  the  querj'  £".V(r^.  z^-  y-,.  y-^  ',;)  is 
handled  in  time  0(log-".  —  .' ).  where  :  is  the  total  number  of  spheres  satiifying  the  query  condi- 
tions. 

This  description  of  the  data  structure  T  completes  the  definition  of  our  algorithm  and  the 
analysis  of  its  complexity.  We  now  turn  to  prove  the  correctness  of  the  algorithm,  h  is  plain  that 
the  algorithm  never  reports  a  spurious  intersection.  Keace  »e  only  need  to  show  that  if  any  two 
spheres  in  S  intersect  each  other,  then  at  least  one  such  intersection  will  be  detected  in  at  least 
one  sweep  performed  by  the  algorithm.  In  fact  we  will  show  thai  the  intersection  between  the 
two  spheres  ^i.  i^  ^°^  which  min(ri.r-;)  is  maximal  mtist  be  detected  in  at  kast  one  sweep  of  the 
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algorithm.  Suppose  that  these  two  spheres  s^,  3o  6  S  intersect  each  other,  and  let  their  centers 
be  Ci,  Co  and  their  radii  be  rj,  r2  respectively,  so  that  ri  <  ro,  and  that  no  two  spheres  both  hav- 
ing radii  greater  than  ri  intersect  one  another  (in  what  follows  we  assume  for  simplicity  that  the 
radii  of  the  spheres  in  5  are  all  distinct).  Let  ^o  €  *  be  a  direction  which  forms  an  angle  6  <  a 
with  the  directed  segment  CjCo,  and  consider  the  sweep  in  the  ?o-direction  which  the  algorithm 
performs.  Rotating  a.xc3  if  necessary,  assume  that  Oq  points  along  the  positive  ;-a.xi3,  and  that  Ci 
is  lower  than  c;.  By  Lemma  1,  the  line  /  parallel  to  the  ;-a.\is  and  passing  through  Cj  intersects 
the  cube  A'  inscribed  in  Sn  with  faces  parallel  to  the  coordinate  planes,  i.e.  the  projection  c/  of  cj 
onto  the  x .y  plane  falls  inside  the  projection  of  K  onto  that  plane.  Hence  it  suffices  to  show  that 
by  the  time  the  swept  plane  P  passes  through  Cj,  the  sphere  Sj  is  still  in  T,  because  then  the  set 
produced  in  response  to  the  range  query  performed  in  step  (a)  of  the  algorithm  at  Co  will  include 
Si.  so  that  the  intersection  between  Si  and  £2  will  be  detected  by  the  algorithm.  Suppose  to  the 
contrary  that  Si  has  been  deleted  from  T  prior  to  the  processing  of  C;.  Then  there  must  exist  a 
third  sphere  s  ^  S  during  whose  processing  Si  was  deleted  from  T.  It  follows  that  the  radius  r  of 
s  is  >  Ti,  that  the  center  c  of  8  lies  above  C;  but  below  co,  and  that  c/  also  lies  in  the  projec- 
tion onto  the  x,y  plane  of  the  cube  A''  associated  with  the  sphere  s.  Moreover,  we  can  also 
assume  that  s  and  s^  do  not  intersect  each  other,  for  otherwise  this  intersection  would  have  been 
detected  by  the  algorithm  at  the  time  c  was  processed.  Also  by  assumption  the  spheres  v'  and  ?; 
do  not  intersect  each  other  because  both  their  radii  are  greater  than  rj.  However,  such  a 
configuration  of  spheres  is  impossible.  Indeed,  let  u  be  a  point  on  /  having  the  same  :  coordinate 
as  Co.    Then  we  have 

c.ii  =  c-ico  sin  e  <  (ri+  ro)  sin  a  <  -'■i^T/j  =  ^ 
Let  J.  J  I.  J2  be  the  intervab  in  which  /  intersects  the  spheres  s,  s^,  ^o  respectively.    Since  .» 
intersects  neither  aj  nor  8^,  and  since  the  center  c  of  ^  lies  above  Cj  but  below  Co,  it  follows  that 
J  is  disjoint  from  both  Ji  and  /o  and  lies  between  these  two  intervab.  It  follows  that 


L«>  \J\  +  Tl-^ii  +  tI-^^I 


But  I  /i]  =  2ri,  I  /ol  =  2>/r2  -  CoU*,  and  |  /|   >  -^  (because  /  intersects  K'  ).  Since  r  >  r, 

v3 


•i"  >  -^  +  '•i  +  v/T7~~^ 


Squaring  the  last  inequality,  we  obtain 


CiU->  rHi+-^)    +    r|-C2u2+   2ral+ -^)V>T^7Mr 
But  CjU*  +    Cou"  =  Cjcl  <  (ri+  r^)^  because  ^j  and  80  intersect  each  other.  Hence 

By  the  first  inequality  above  we  have 

>yri  -  c,u^  >  s/7rT77W=  ^- 

Substituting  this  in  the  preceding  inequality,  and  dividing  by  rj,  we  obtain 


r,  +  2r2> 


li^^l-h^l^'- 


which  is  impossible,  since  the  coefficients  of  Tj  and  of  rj  are  both  positive.  This  contradiction 
implies  that  the  assumed  configuration  of  spheres  is  impossible,  which  in  turn  implies  that  the 
intersection  betv/een  ?i  and  So  will  be  detected,  proving  the  correctness  of  our  algorithm,  thus  jus- 
tifying the  following  summary  theorem. 

Theorem:  The  existence  of  an  intersection  among  n  arbitrary  spheres  in  S-space  can  be  tested 
by  an  algorithm  requiring  0(n  log  n)  space  and  0(n  log^n)  time. 

Remark:  The  algorithm  described  above  can  easily  be  generalized  to  any  number  d  >  3  of 
dimensions.  The  appropriate  generalization  of  the  data  structure  T  is  a  (<f-2)-dimensional  range 


.'3r;t:ri(Tr'0  c.   tTTi'df'-!'- 
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tree  with  a  priority  search  tree  attached  to  each  of  its  (innermost)  nodes.  The  set  *  of  directions 
in  which  sweeps  are  to  be  performed  must  be  dense  enough  so  that  for  each  orientation  0  in  the 
((/-l)-dimen3ionai  unit  sphere  S'''^  there  exists  an  orientation  0  6  *  such  that  the  planar  angle 
between  6  and  ip  is  <a^,  where  sin  a^  =  — —.  The  correctness  proof  uses  straightforward  gen- 
eralizations  of  the  calculations  given  above,  in  which  VS  is  replaced  by  \/J.    This  gives 

Corollary:  The  existence  of  an  intersection  among  n  arbitrary  rf-spheres  in  Euclidean  rf-space 
can  be  tested  by  an  algorithm  requiring  0(n  log''"^n)  space  and  0{n  log''"'n)  time. 
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